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Abstract 

First principles calculations based on the density functional theory within the local 
spin density approximation plus U(LSDA+U)scheme, show rhombohedral Bi2FeTi06 
is a potential multiferroic in which the magnetism and ferroelectricity coexist . A 
ferromagnetic configuration with magnetic moment of 4 /x^ per formula unit have 
been reported with respect to the minimum total energy. Spontaneous polarization 
of 27.3 /i C/cm^ ,caused mainly by the ferroelectric distortions of Ti, was eval- 
uated using the berry phase approach in the modern theory of polarization. The 
Bi-6s stereochemical activity of long-pair and the 'd^-ness' criterion in off-centring 
of Ti were coexisting in the predicted new system. In view of the oxidation state of 
Bi^+ jFe^"^ ,Ti^+ , and from the orbital-resolved density of states of the Bi-6p, 
Fe-3d,Ti-3d , and 0-2p states, the valence state of Bi2FeTi06 in the rhombohedral 
phase was found to be Bia^+Fe^+Ti^+Oe. 
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1 Introduction 



Multiferroic materials liave coupled magnetic, electric, and even structural or- 
der parameters that result in ferromagnetic, ferroelectric, and/or ferroelastic 
properties at same phase. The magnetization can be rotated or even reversed 
by the ferroelectic reversal, and vise versa. These materials have potential ap- 
plications in information storage, the emerging field of spintronics,and sensors. 
Since 1960s various articles have systematically classified their properties and 
behaviors [T|2|3llHl5] . and the theoretical study was also carried out on the mag- 
netoelectric coupling effect [6. 7J. BiFeOs was the rare one in nature, which pos- 
sess both weak ferromagnetism and ferroelectric characteristics in single phase. 
Meanwhile, more experimental investigation have been performed on the mag- 
netoelectric properties of it . BiFeOa has long been known to be ferroelectric 
with a Curie temperature of about 1103 K and antiferromagnetic(AFM) with 
a Neel temperature of 643 K. The Fe magnetic moments were coupled fer- 
romagnetically in the (111) plane and antiferromagnetically in the adjacent 
plane along [111] crystalline direction, which is called the G-type AFM. The 
rhombohedral distorted perovskite structure with space group R3c permits 
a canting of AFM sublattice, resulting in a weak ferromagnetism. However, 
there is a spiral spin structure in which the antiferromagnetic axis rotates 
through the crystal with a long- wavelength period of 620A [8J. Therefore, the 
linear magnetoelectric coupling was destroyed with the cancellation of mag- 
netization, due to this long spiral spin structure. Although this cancellation 
could be suppressed partly in thin film , the magnetization observed in the 
BiFeOa thin film was relatively smaller in terms of the excellent ferroelectric 
behavior. 
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This cancellation of magnetization might be inhibited by doping foreign mag- 
netic transitional metal ions such as Cr,Mn,Co,and Ni in the B sites. BiCrOs 
was recently synthesized and reported to be a highly distorted perovskite with 
C2 symmetry [9j, which permit the occurrence of ferroelectricity. Unfortunately 
, antiferromagnetic ordering was reported with a residual magnetization of 
0.02 /i^jSuggestive of weak ferromagnetism. Bi2FeCr06 has been predicted 
theoretically to be ferrimagnetic(2 per formula unit)and ferroelectric(^80 
/iC/cm^). The predicted ground state structure is very similar to the _R3c struc- 
ture of BiFeOa, except that the Fe^"*" in the second (111) plane are replaced by 
the Cr^"*", which lead to the space group i?3. 

It is known that the ferroelectricity in BiFeOs is caused by the long Bi-6s 
stereochemically active pair, which formed the one chemical mechanism for 
stabilizing the ferroelectricity . This mechanism was caused by the mixing be- 
tween the (ns)^ ground state and a low-lying {nsYinpY excited state, which 
can only occur if the cation ionic site does not have inversion symmetry. The 
second mechanism leading to the polarization in ferroelectrics is the off-centre 
displacement of small cation in the common perovskite ferroelectrics such as 
BaTiOs, which is attributed to the ligand- field hybridization of the transi- 
tional metal cations by its surrounding anions. This mechanism requires the 
'd^-ness' criterion in the small transitional metal off-centring, hence precludes 
the coexistence of magnetism and ferroelectricity. Sometimes these two mech- 
anisms occur in the same material ,such as PbTiOs. 

Due to the theories mentioned above, the Ti-doping BiFeOs could maintain 
these two behaviors together in same phase and generate the coexistence of 
magnetism and ferroelectricity. The new modified material Bi2FeTi06 with 
rhombohedral structure was expected to have excellent magnetic and ferro- 
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electric properties. As in Bi2FeCr06, the Fe^"*" in the second (111) plane are 
replaced by the Ti'^^, which lead to the space group R3. Up to now, there 
is no theoretical prediction on this new system. In this paper we have sys- 
tematically investigated the electronic properties, ferroelectric behaviors, and 
magnetic ordering of the new theoretically predicted Bi2FeTi06 , using first 
principles calculation based on the density functional theory. 

The remainder of this article is structured as follows:In section II, we presented 
the computational details of our calculations. We provided the calculated re- 
sults and discussions in section III. In section IV , the conclusion based on our 
calculation were given. 

2 Computational details 

2. 1 Methodology 

Calculations in this work have been done using the Quantum-ESPRESSO 
package [TU], which is based on the density functional theory (DFT) and density 
functional perturbation theory(DFPT) using the plane-wave pseudopotential 
formalism. We used our self-interaction-corrected ultrasoft pseudopotential 
implementation with the Perdew Burke Ernzerhof (PBE) exchange correla- 
tion functional, as the common local density approximation fail to obtain a 
band gap in the transitional metal oxides. In our DFT computations, plane- 
wave basis set with kinetic energy cut-off of 40 Ry was employed, and it has 
shown that the results are well converged at this cutoff. Bi 5d ,6s, and 6p 
electrons, Fe 3s, 3p, and 3d electrons, Ti 3s, 3p, and 3d electrons, and O 2s and 
O 2p electrons have been treated as valence states. We used up to 6 x 6 x 6 
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grids of special k points in total energy calculation. In the local spin density 
approximation plus U(LSDA+U) framework, the strong Coulomb repulsion is 
treated by adding a Hubbard-like term to the effective potential, leading to an 
improved description of the interaction in the transitional metal oxides. The 
Hubbard parameter U in the range eV - 8 eV were used in the calculation 
of electronic structures. Results show that the U of 7 eV is enough to produce 
a better insulated solution, which is used in the following sections. 

2.2 Structure optimization 

The lattice parameters was evaluated with respect to the lowest total energy 
, ranging from 5.409 to 5.549 A. In this process the structure was constrained 
to be rhombohedral within i?3 space group. The atomic positions were taken 
from the corresponding value of BiFeOa in Ref . [llj , as the Fe^"*" in the second 
(111) plane are replaced by the Ti^"*" . The calculated lattice constant value 
was plotted with respect to energy and fitted to the Murnaghan equation of 
states ,which give the bulk moduluspjj. The atomic positions were relaxed af- 
terward by minimization of the Hellman-Feynman forces within a convergency 
threshold of 10~^ Ry/Bohr , using Broyden-Fletcher-Goldfarb-Shanno(BFGS) 
variable metric method, as cell shape and volume was fixed. 

We doubled the size of the unit cell along one of the rhombohedral lattice 
vectors and compared the total energies for all possible coUinear spin configu- 
rations. The ferromagnetic (FM) and antiferromagnetic(AFM) structure (with 
nearest neighbor Fe antiparallel to each other) were constructed within the 
resulting supercell up to 20 atoms. Calculated energy show that the AFM 
structure supercell possess a relatively high energy of 0.03 eV per unit cell in 
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comparison with the FM one. The FM ordering along z and x(y) direction 
were considered to find the stable FM configuration. Total energy per formula 
unit show the FM ordering along z gave the stable magnetic configuration. 
Hence, in the following computation we took the FM ordering along z as the 
stable state . 

2.3 Spontaneous polarization 

In the modern theory of polarization approach [T311 14| 1 5j . the total polarization 
P for given crystalline geometry can be calculated as the sum of ionic and 
electronic contributions. The ionic contribution is calculated by summing the 
product of the position of each ion in the unit cell with the nominal charge 
of its rigid core. The electronic contribution is determined by evaluating the 
phase of the product of overlaps between cell-periodic Bloch functions along 
a densely-sampled string of neighboring points in k-space. 

3 Results and discussion 

The optimization of volume was performed with lattice constant ranging from 
5.409 to 5.549 A, and the volume plotted versus energy was shown in Fig.l. 
The calculated value was fitted to the Murnaghan equation of states, per- 
mitting to evaluate the bulk modulus. The normal parabolic curve gave the 
optimal value of 5.506A , compared with the corresponding value of 5.459A for 
BiFeOa^lI]. We also reported the magnetic moment per formula unit versus 
volume so as to see the variation of it with respect to volume change, and the 
optimal volume give a magnetic moment of 3.99 per formula unit. Next, 
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the atomic positions were relaxed using the optimal lattice constant with cell 
shape and volume fixed. The relaxing Wyckoff positions, calculated lattice 
constant ,rhombohedral angles, and unit cell volumes were all summarized in 
Table I. The internal structural parameters are very similar to those in BiFeOa 
,except that the O positions were deviated from the initial positions. This 
may partly reflects the fact that there were change of hybridization in the new 
predicted BiaFeTiOe- 

For the simple system consisting of N atoms, transforming the total Humilto- 
nian into Wannier representation for the narrow band gap condition we get 

H = Y^Y^Tr^CfsC^s + ^ E E < ^Mlm > C+C+^CmS'Qs (1) 

i,j 5 i,j,l,m S,S' 

The second term that is the interaction part can be simplified with considering 
the short distance effect on one lattice or atom as: 

E E c+c+Qs'Qs = luY.E (2) 

^ i S,5' ^ i S 

riis and n^g is the operator of spin-up and spin-down number of particles. 

Here U is 



2 /■«*(r-Ri)a*(r'-Ri)a(r-Ri)a(r'-Ri) , 

<ii\v\ii>—e / — ^ — ^ 1 — -drdr (3) 

J |i" ~ I" I 

Hence, U is the Coulomb interactions between the electrons around the lattice 
or atom and take great effect on the metal-insulator transition. The electronic 
structures were calculated with U varying from eV to 8 eV. The results 
show the U of 7 eV is able to give a better insulated solution, while the 

insulated characteristics have not been improved significantly after increasing 
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of U further. We used this value of U in the following sections. The total 
density of states(DOS) for Bi2FcTi06 before and after applying U was shown 
in Fig. 2. It can be seen the band gap was increased from 0.46 eV (U=0 eV)to 
1.68 eV(U=7 eV). The top valence bands caused by Ti was pushed to the 

conduction bands after employed the U of 7 eV. 

In order to gain more insight into the role of ferroelectric distortions on elec- 
tronic structure and bonding behavior in these new predicted Bi2FeTi06. The 
calculated band structure of up-spin states was given in Fig. 3. There are strong 
hybridization between Bi-6p and 0-2p states in the top of the valence bands, 
and the Fe-3d states was pushed down in lower energy region compared with 
those in BiFeOa. The strong hybridization also reflect the fact that the activity 
of Bi-6s long pair also take place in this new R3 structure. The band at -10 eV 
are mainly originating from the Bi-6s electrons ,these two bands are relative 
broader in this new system. Ti-3d bands mainly positioned at the bottom of 
conduction bands and hybridized deeply with 0-2p states in this region, and 
the Ti 3d also overlap shghtly with O 2p under the Fermi energy. 

To understand the change of hybridization of Ti-0 and Fe-0 after Ti doping 
, the calculated charge density in the energy window from Ep-G eV to Ep 
covering the region ,where the Ti-0 bonding states are located ,were shown 
in Fig.4 and Fig.5,respectively. From these two Figures, it is clear that Ti is 
bonding strongly with O ,while Fe-0 bonding strength is reduced comparing 
with those in BiFeOa. Moreover the Ti-0 bonding is along the direction in 
which the Ti-3d orbitals formed and slightly affected by the hgand-fleld. It 
means the Ti is bonding deeply with the nearest six Oxygen atoms. Therefore 
the Ti off-centring driving ferroelectric displacement still maintained in these 
new system. We are expecting this chemical mechanism for stabilizing the 
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distorted structure in ferroelectric oxides take effect in tfie new system. 

Tfie Ti tend to form the oxidation state of Ti^+, which can be approved by 
the orbital-resolved DOS for Ti shown in Fig. 6 consequently. It is found the 
dxy, dyz, , dxz,Sirid dx2-y2 orbitals , including the majority spins and minor- 
ity spins, are all nearly unoccupied in the conduction bands indicating the 
oxidation state of Ti'^^. These bands arc hybridizing strongly with the 0-2p 
electrons, leading to the stabilizing ferroelectric off-centring behavior. And this 
phenonmenon can further confirm the occurrence of this chemical mechanism 
in this predicted new system. The dyz, dxz,sind d^^-y-i orbitals have rather little 
states broadened under the Fermi energy, which indicate the smaller deviation 
of oxidation state from Ti^+. These bands also overlap with 0-2p electrons to 
strengthen the hybridization. 

In the new structure, Fe atom tend to form oxidation state of Fe^"*" ,which is 
lowered comparing with the Fe^+ in BiFeOa. We reported the orbital- resolved 
DOS for Fc in BiFeOs and Bi2FcTi06 shown in Fig. 7 and Fig.8 respectively , so 
as to observe the change of oxidation states. It can be seen in Fig. 7 that all the 
majority spin states are occupied while all the minority spin states are nearly 
unoccupied with fewer states lying at -3.5 eV , indicating the oxidation state 
of Fe^+ in BiFeOs. However ,from orbital-resolved DOS of Fe in Bi2FeTi06 
shown in Fig.8 , unlike the case in BiFeOs ; the spin-down states in (lj.y and 
dxz orbitals are partly occupied , indicating the orientation of bonding in this 
circumstance. Meanwhile all the spin-up states are occupied as in BiFeOs. 
Based on the large amount of occupied spin-down states in d^y orbitals and 
the relatively fewer occupied spin-down states in dxz orbitals, we deduced the 
formal spin (S=4/2) of Fe^+ with d^ configuration in Bi2FeTi06. Although the 
Fe-3d states were narrowed and pushed close to the Fermi energy after doping 
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with Ti, the main states have not overlapped deeply with the 0-2p electrons. 
Hence the Fe 3d-0 2p hybridization have not been strengthened consequently 

The Bi-6s stereochemically active long pair were participating in the ferroelec- 
tric distortions. The hybridization between Bi-6p electrons and 0-2p electrons 
can be observed in the orbital-resolved DOS of Bi and O in Bi2FeTi06 shown 
in Fig.9. We also reported the total DOS in Fig. 9 for comparision. It is clear 
that the Bi-6p states and 0-2p states positioned at the same energy region 
and hybridized strongly as in BiFeOa. Moreover,the Bi 6p does not change 
much comparing with the states in BiFeOa , indicating the similar oxidation 
states in these two systems. In view of the oxidation states of Ti^+,Fe^+,and 
Bi3+, the valence state in the R3 BiaFeTiOg is Big^+Fe^+Ti'^+Oe. 

The spontaneous polarization was evaluated using the berry phase method 
based on modern theory of polarization. And Born effective charges(BECs) 
were estimated by computing the Cartesian components of the polarization 
with respect to the atomic displacements, ie.: 




where A//j^ is the displacement of ion j in Cartesian direction /3, Z^^^ is its 
BECs tensor, and VL is the unit cell volume. The P is calculated from berry 
phase method by considering a specific structural pathway parameterized by 
the change in atomic displacement connecting a centrosymmetric reference 
structure and i?3. The displacement along [111] direction was chosen smaller 
(2% lattice constant)enough to ensure the vahdity of the hnear treatment in 
Eq.(4). The results were listed in table II. 
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The new predicted Bi2FeTi06 was found to obtain a medium value of spon- 
taneous polarization of 27.3 fi C/cm^ , comparing with the relative large value 
of 59.4 fi C/cm^ in BiFeOa. This result was not so good as we expected in 
terms of the coexistence of these two chemical mechanisms of ferroelectric 
distortions in the same phase. This phenomenon might be partly explained 
by the anomalous large EEC of 5.82 for Ti , leading to the dominant role of 
Ti off-centring distortions in the new predicted system. The producing polar- 
ization was similar to the value in BaTiOs [16j in which the same driving 
mechanism applied. Moreover, the EEC of Bi and Fe were all decreased in the 
new system, indicating the decreasing ability of producing the ferroelectric- 
ity of Bi-6s long-pair. There might be interaction between these two chemical 
mechanisms for producing ferroelectricity , which need to be explored further. 
The predominant role of Ti off-centring ferroelectric displacement generated 
a smaller ferroelectricity than the Bi-6s long-pair in terms of a collapse of 
32 yU C/cm^ in Bi2FeTi06 compared with that in BiFeOs . Fortunately, the 
adoption of Ti in the new system inhibited the cancelling of macrosopic weak 
magnetization , accompanying the oxidation state of Fe^"*" with configu- 
ration(with magnetic moment of 4 fiB per formula unit ). Consequently, the 
coexistence of ferromagnetism and ferroelectricity can be fulfilled in the new 
predicted Bi2FeTi06. 

4 Conclusion 

The cancellation of magnetization caused by the spiral spin structure order- 
ing in BiFeOa can be improved by doping foreign transitional metal ions. 
Bi2FeTi06 was predicted to be a potential multiferroic in which a large fer- 
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romagnetism(with magnetic moment of 4 //b per unit cell) and ferroelectric- 
ity coexisted. Berry phase calculation show a medium polarization of 27.3 /j, 
C/cm^ in Bi2FeTi06 comparing with the value of 59.4 n C/cm^ in BiFeOs, 
which was mainly caused by the Ti off-centre ferroelectric distortions, although 
the hybridization between Bi-6p states and 0-2p states had taken place in this 
new predicted system , which was reduced after doping Ti , accompanying a 
small deviation of BECs of Bi from the nominal one. Prom the orbital-resolved 
DOS of the elements in the system ,the oxidation state of Bi2FeTi06 was found 
to be Bi|"'"Fe^"'"Ti^+06 with decreasing of oxidation state of Fe from Fe^+ to 
Fe2+. 

This work was supported by the Aeronautical Science foundation of China 
(Grant Nos. 2003ZG51069). 
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Table 1 

Calculated lattice constant a, rhombohedral angle a, volume V, and Wyckoff 
parameters for Bi2FeTi06. The WyckofF positions la(x,x,x) for the cations and 
3b(x,y,z) for the anions. 

Bi2FeTi06 BiFeOs 

5.506 5.459 

60.36 60.36 

118.97 115.98 

X -0.001/0.501 0.000 

X 0.229(Fe)/0.728(Ti) 0.231 

X 0.558 0.542 

y 0.935 0.943 

z 0.393 0.398 

Table 2 

The polarization (P) and Born effective charges(BECs)(Z*) for displacements along 
[111] for Bi2FeTiQ6 and BiFeOs . 
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Figure captions: 

Fig.l Volume vs. total energy per formula unit and magnetic moment per 
formula unit for Bi2FeTi06 . 

Fig.2 (a) and (b) show total DOS for BiaFeTiOe with U=0 eV and U=7 eV 
respectively. The Fermi level was set to zero. 

Fig. 3 Electronic band dispersion for Bi2FeTi06. The Fermi level was set to 
zero. 

Fig. 4 Charge density (in arbitrary units) in a particular plane for BiFeOs. 
The charge density is calculated in the energy window from Ep-G eV to Ep. 

Fig. 5 Charge density (in arbitrary units) in a particular plane for Bi2FeTi06. 
The charge density is calculated in the energy window from Ep-Q eV to Ep. 

Fig. 6 Orbital-resolved DOS for Ti in Bi2FeTi06.(a),(b),(c),(d),and (e) show 

the DOS for d^y, dyz, dz-2 , dxz,s^-nd d^^^y'z orbitals rcspcctivcly.Spin-up states 
are shown in the upper portions and spin-down states in the lower portions 
in all panels.The Fermi level was set to zero. 

Fig. 7 Orbital-resolved DOS for Fe in BiFe03.(a),(b),(c),(d),and (e) show the 
DOS for dxy,dyz,dz2,dxz, and d^^-y^ orbitals respectively. (f) show the Fe 3d 
total DOS. Spin-up states are shown in the upper portions and spin-down 
states in the lower portions in all panels.The Fermi level was set to zero. 

Fig. 8 Orbital-resolved DOS for Fe in Bi2FeTi06.(a),(b),(c),(d),and (e) show 
the dxy,dyz,dz2,dxz,Q'ri(idx2_y2 orbitals respectively. Spin- up states are shown 
in the upper portions and spin-down states in the lower portions in all 
panels.The Fermi level was set to zero. 
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Fig. 9 Total DOS and orbital-resolved DOS for Bi 6p and O 2p in 
Bi2FeTi06.(a) and (b)sliow the Bi-6p and 0-2s states respectively, (c) is the 
total DOS. Spin-up states are shown in the upper portions and spin-down 
states in the lower portions in (c).The Fermi level was set to zero. 
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This figure "Figure4.jpg" is available in "jpg" format from: 
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This figure "Figure5.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/0710.2391vl 
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